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Abstract 

We propose a new scheme for parallel spatially multimode quantum memory for light. The 
scheme is based on the propagating in different directions quantum signal wave and strong classical 
reference wave, like in a classical volume hologram and the previously proposed quantum volume 
hologram [IJ. The medium for the hologram consists of a spatially extended ensemble of cold 
spin-polarized atoms. In absence of the collective spin rotation during the interaction, two passes 
of light for both storage and retrieval are required, and therefore the present scheme can be called 
a double pass quantum volume hologram. The scheme is less sensitive to diffraction and therefore 
is capable of achieving higher density of storage of spatial modes as compared to the thin quantum 
hologram of [2], which also requires two passes of light for both storage and retrieval. On the 
other hand, the present scheme allows to achieve a good memory performance with a lower optical 
depth of the atomic sample as compared to the quantum volume hologram. A quantum hologram 
capable of storing entangled images can become an important ingredient in quantum information 
processing and quantum imaging. 
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I. INTRODUCTION 



A number of quantum information protocols such as quantum repeaters, distributed 
quantum computation, quantum networks etc. require or will greatly benefit from using a 
quantum memory. A variety of approaches for storage in atomic ensembles were developed 
recently, including the schemes based on quantum nondemolition (QND) interaction, elec- 
tromagnetically induced transparency (EIT), the Raman scattering, and photon echo. A 
comprehensive recent review on quantum interfaces between light and matter can be found 
in [3]. The problem of multimode quantum memories is at the center of current research due 
to their potential for enhanced storage capacity, which is necessary for scalable linear-optical 
quantum computing [4J and efficient quantum repeaters [H]. 

In this paper we present a new scheme for the spatially multimode quantum memory 
for light, which we call a double pass quantum volume hologram. Our proposal, like the 
previously considered quantum volume hologram [Ij, makes use of the concept of counter- 
propagating geometry which stems from a classical volume hologram, one of the corner 
stones of modern holography |6]. When a hologram is written by the counter-propagating 
or, in general case, propagating in different directions signal and reference waves, the two 
sublattices produced by the waves interfering in the medium are recorded, each of them 
storing one quadrature of the signal field. Since both quadratures are stored, there are no 
virtual and real images during the readout, in contrast to a classical hologram recorded in 
a single pass of the co-propagating waves. 

In the counter-propagating geometry there is no need for phase matching between the 
signal and the strong reference waves because their relative phase oscillates in space hundreds 
or thousands times. This feature significantly weakens [1] the diffraction limitation on the 
transverse spatial density of stored modes as compared to the co-propagating geometry of 
the thin quantum hologram of [2] . 

A limitation of the quantum volume hologram with rotating spins llj and of the closely 
related Raman memories for light, extensively studied in [7HTT]. is due to their single-pass 
operation. The state exchange between light and atoms within the propagation length 
is associated with a "self-erasing" of quantum state of the input light, which is, roughly 
speaking, exponential with oscillations in the longitudinal direction. The higher efficiency 
of the state exchange is needed, the larger optical depth of the sample should be provided. 
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In contrast to [T], we consider here the ground state atomic spins in absence of rotation 
during the interaction. Hence, during one pass of hght the longitudinal (with respect to the 
signal wave propagation direction) quadrature component of the collective spin in a given 
sublattice does not change its state and, at the same time, is effectively recorded into one 
of the signal wave quadratures within all duration of the interaction cycle. This feature of 
our present scheme is inherited from the quantum non- demolition interaction, typical for 
the thin quantum hologram, and allows for a good efficiency of the write-in and read-out by 
a fixed value of the coupling constant, which is smaller than for the the previously proposed 
quantum volume hologram or the Raman memories. 

Due to the non-demolition in part character of the present scheme, two passes of light are 
needed both at the write-in and the read-out stages in order to "erase" initial quantum state 
of the light signal and the atomic spins by the state exchange, similar to the thin quantum 
hologram. 

We evaluate the transverse spatial density of the stored light modes, and calculate the 
average quantum fidelity per transverse mode (per pixel) for the whole storage cycle. For 
the initial vacuum (non-squeezed) state of collective spin, the upper limit on the fidelity is 
given by Fav ~ 0.845, and the fidelity can be made close to 1 given an effective initial (before 
the write-in stage) squeezing of first longitudinal modes of the collective atomic spin. 

II. SINGLE PASS OPERATION OF VOLUME HOLOGRAM 

Consider an ensemble of motionless atoms with spin 1/2 both in the ground and in the 
excited state, located at random positions. The long-lived ground state spin J° of an atom 
is initially oriented in the vertical direction x. A classical off-resonant x-polarized plane 
wave at frequency ujq with a slowly varying amplitude (assumed to be real) propagates 
in the direction fcc, where the corresponding wave vector kc lies in the (y, z) plane. In 
this paper we investigate the evolution of an input signal which is represented by a weak 
quantized y-polarized field at the same frequency Wq; propagating in +z direction. In what 
follows we consider this spatially multimode input field with a slowly varying amplitude 
Ay{f,t) ^ Ax in the paraxial approximation. In our scheme the resulting field strength 
oscillates in the (x, y) plane. Let us introduce the difference wave vector Ak = kg — kc 
with the only non-zero projections Aky and Ak^- One can consider a thin (that is, of the 
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Hologram 




FIG. 1. Illumination of the hologram by the write-in (single pass operation). 



width much less than A) atomic layer orthogonal to Ak, where the phase difference between 
the signal and the driving classical wave is constant. The quantum non-demolition (QND) 
light-matter interaction in each layer leads to two basic effects: (i) the Faraday rotation of 
light polarization due to quantum 2;-component of collective atomic spin of the slice; and (ii) 
the atomic spin rotation, caused by the unequal light shifts of the ground state sub-levels 
with rriz = ±1/2 in the presence of quantum fluctuations of circular light polarizations. The 
interaction within the slice is described by the well known QND Hamiltonian [3j: 

2nko\d\^ 



H 



eg 



I dfY^j:{t)Sz{f,t)6{f-f,) 



(1) 



Here Ueg is the frequency of the atomic transition, d is the dipole matrix element, and 
ko = uqc. For the non-copropagating signal and driving waves, the z-component of the 



Stokes vector Sz{r,t) 



2A^ Im 



Ay{f,t) exp{iAkr) in (1) slowly varies within the slice 
and rapidly oscillates (i.e. on the scale 1/Ak) along the Ak direction. The slowly varying 
amplitude Ay of the quantized signal field is defined via 

dky, f dq 



Ay{z,p,t) 



uj{k)/ko ay{k) exp[i{qp+ {k^ - ko)z - {u{k) - uo)t}], 



27r J (27r)2 

where f = {p,z), k = {q,kz). Here ay{k) and a^y{k) are the annihilation and creation 
operators for the wave k, which obey standard commutation relations [ay{k) , a'^y{k ')] = 
(27c)^6{k — k'), [ay{k) , ay{k ')] = 0. By using these commutation relations in the paraxial 
approximation, one finds [I2] the commutation relation for the slowly varying amplitude of 
the quantized signal field Ay, 



[Ay{z, p, t), AUz', p', t)] = c5(f - f ') 



(2) 



We do not consider here quantized waves of other than kg mean propagation directions 
because their evolution is independent of the signal wave under consideration. 

We introduce the density of the collective spin as J(r) = J2a J"'5{f — Ta). The averaged 
over random positions of the atoms commutation relation for the y, z components of the 
collective spin is 



r — r") 



Here is the average density of atoms. The field-like canonical variables for the spin 

subsystem, 

X{r,t) = Jy{f,t)/^na{JS), P(r,t) = J,(f, t)/vt(J^, 
obey the canonical commutation relation: 



[X{f, t), P{f', t)] = i6{f- f '). (3) 

The full Hamiltonian of our model includes the energy of free electromagnetic field and the 
effective Hamiltonian of QND interaction. The Hamiltonian reads, 

H = J^df!^^Al{z,p,t)Ay{z,p,t)- (4) 

.27rko\d\ 



-./^)A,{z, t)P{z, p, t) [Ay{z, p, t)e'^^''^ - h.c] ] . 



We describe the evolution of our system in the Heisenberg picture. With the use of commu- 
tation relations ^ and ^ for the field and atomic variables, after simple transformations 
we obtain: 



<9 ^ . 2R 



9t ^ ^ yzr 
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X{z, p, t) = -j^lm Ay{z, p, t)e'^'^' , -P{z, p, t) = 0. (6) 



Here T is the duration of the flat-top pulse of driving field and L is the atomic cell length. 
The dimensionless coupling constant 

^= ,';^°''[, v^.(^gM^^T, (7) 

should be of the order of unity for the memory to work (note that R is a/2 times the coupling 
constant k defined in [1] for volume hologram with rotating spins). The coupling constant 



can be written as = 2ao^, where ao is the resonant optical depth and t] is the probabihty 
of spontaneous emission [13]. Since 1] 1 is required in order to neglect the effect of 
spontaneous emission from the upper level, the usual condition = X'^naL/2n ^ 1 should 
be fulfilled. We introduce the Fourier transform via 



a{z,q,t)= dpAy{z,p,t)e~ 



■iqp 



(8) 



and similar for atomic variables, and arrive to the set of basic equations in the Fourier 
domain: 



d 



dt 
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dz 2kQ 
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p{z,q- Aky,t)e~ 



iAk,z 



2k 



Jm 



d 



dt 



p{z,q,t) = 0. 



(9) 
(10) 



a{z,q + Aky,t)e'^''^'\ , 

Notice that the collective spin quadrature amplitude p responsible for the Faraday rotation 
does not evolve, similar to the case of thin quantum hologram of [2j. Hence, this amplitude 
is effectively recorded into one of the signal wave quadratures within all duration T of one 
interaction cycle. 

Let us consider the following boundary condition. We position the center of atomic 
sample of length L at z = 0, and define the input signal wavefront A*^*"^(p, t) and the Fourier 
amplitudes a''^'^\q,t) with respect to the central (x, y) plane at z = in free space. The 



actual signal field at the input cell face is related to the defined "in" amplitudes as 

a(-L/2,g,t) = a^''^\q,t)ei^ . 



(11) 



One can imagine an external lens focusing system which transfers the field A*^*"^ (p, t) from 
its input plane to the plane at 2; = 0, so that A*^*")(p, t) is the input field of the lens system. 
We arrive to the following solution of the Eq. ([9]): 



K 



dz'p^'''\z',q- Aky)e- 



-iAk^z' 



-iz-z') 



(12) 
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where p^^"'\z,q) = p{z,q,0). We assume pulse duration T to be much larger than the 
retardation time at the length L of ensemble, T ^ L/c and neglect ~ 1/c term in the field 
evolution equation The Eq. (10) yields, 

2k 



x{z,q,t) = x^''\z,q) + 



dt' Im 



a{z, q + Aky, t )e 



l\ AAk, 



(13) 
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We should have in mind that the field amplitudes are defined as slow in space in dependence 



of z coordinate, but the spin amplitudes are not. The Eqs. (12), (13) show that QND 



interaction with quantized field a(z, q^t) produces fast spatial modulation of collective spin, 
and similar fast modulation is read out by the the field. The fast modulation of collective spin 
at spatial frequency A/c is just a consequence of the non-collinear geometry of our volume 
hologram. The volume hologram is a spatial multi-layer structure with typical spatial period 
of the order of 27r/Afc. Within this spatial period, the phase difference between the signal 
and the driving classical wave changes by 27r and therefore changes the type of local circular 
polarization which is crucial for QND interaction. This imposes limitations on the atomic 
motion during storage time: the atoms should not transport coherence from one layer to 
another. A solid state, an optical lattice, or ultra cold atoms should be used. 

We go over to slow amplitudes of the collective spin by compensating fast spatial oscilla- 



tions present at right side of Eqs. (12), (13), and perform mode decomposition of the slow 
amplitudes: 



rL/2 2 

/-L/2 

The same decomposition is applied to quadrature amplitude p{z, q — Aky). Here 



xMt)= I dze^{z)xiz,q- Aky,t)e-'^'^''^~^o^\ (14) 

J-L/2 



^n(^)=Ar„^|p.(|), ^n = f-^, (15) 

Pn{x) is the n-th Legendre polynomial, n = 0, 1, 2, . . ., and Afn is the normalization constant. 
A similar mode decomposition in time domain has been introduced in Some lower-order 
amplitudes are given by 

~1 rL/2 2 

L J-L/2 



= \Iy r" dzxiz,q~Aky,t)e '^^"^ , (16) 



xMt) = 2 ( -^j J^^^^dzzx{z,q- Aky,t)e-'''^''-^o'\ (17) 



x,{q,t) = 6 {^^^'^ f_[\dz ^i^,q- Aky,t)e-'^^'^-i^\ (18) 
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We assume that the length L of atomic cell is large as compared to the wavelength. For 
large enough angle between the driving and the signal field wave vectors kc and kg (that is, 
a collinear geometry is not considered here) one has Ak^L ^ 1. Physically this means that 
our volume hologram contains many interference layers discussed above. 

These new amplitudes allow us to express solutions for atomic and light variables in a 
more compact way. As before, we define the output signal field amplitude with respect to 
the central {x,y) plane at 2; = 0. The actual signal field at the output face of the cell is 

a(L/2, g, t) = a^""*) (g, t)e~'i^ . (19) 



Here one can imagine an output lens focusing system which in free space transfers the signal 
field from the plane z = to its output plane. Define the averaged over the pulse duration 
T signal amplitude as 



a{q[ 



Vt Jt 



The field propagation equation ( 12 ) yields 



a(°«*)(g-') = a^-^(g') + ^pr(g). 



(in) 



dt a{q, t). 



(20) 



(21) 



In order to find slow atomic amplitudes we insert ( 12 ) into ( 13 ) and make use of the definition 



(14). The proportional to e '^^^^^ term in (13) after substitution into (14) gives rise to fast 



spatial oscillations, and in the limit Ak^L ^ 1 vanishes by the integration over z. For the 
output slow amplitudes of the collective spin, where x^°"*^(g) = Xniq,T) and -*((?) = 
Pniq,T), we find 

1 



xr'\q) = xr\q) - ika^'''\q) - i 



(in) / ->\ J- (in) / ->\ 

~ 71 



where tt, = 1, 2, 



and 



(in). 



Here the matrix Q is given by 



Qn 



71-1 



(2n- l)(2n + 1) 



Qn,n 



+1 



(2n + l)(2n + 3) 



(22) 
(23) 

(24) 
(25) 



While deriving (22, 23j) we made use of the following property of the Legendre polynomials, 

J dxPn{x) = ^^^[Pn+l{x) - Pn-l{x)]. (26) 



The Eqs. (21, 22) demonstrate that by a single pass of light, the amplitude Po"'\q) of the 
collective spin is read-out by the signal, and the amplitude a^™)(g) of signal field is mapped 



onto the atomic amplitude Xo°"*''((f). One can also observe in the Eq. (22) effects of the 
transfer of the collective spin coherence between atomic layers in the positive 2;-direction, 
which are absent in the QND based thin hologram. The "active" spin amplitude p produces 
a contribution to the quantized field a, which is written into the "passive" amplitude x 
somewhere in the next layers. 
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By contrast to the QND based scheme of p], both quadrature amphtudes of the hght 
field are stored in a single cycle of light-matter interaction. This feature of the present 
scheme is related to its non-collinear geometry. One can imagine two shifted by a fraction 
of wavelength sub-lattices of the collective spin, so that any of two quadrature amplitudes 
of the signal is independently stored in the relevant sublattice. 

In some sense a single-pass storage and a single-pass retrieval of the signal light state is 
achieved in the present scheme. After the first (the write-in) cycle of interaction one has to 
rotate the spins by 7r/2 relative the x-axis, so that the initial condition for the second (the 
read-out) cycle is given by 

p^^'^'Hq) = x^^'^-'\q), (27) 
and let another light pulse to read-out the stored state. By applying twice the in-out 



transformations (21, 22), we find for the output light amplitude: 

7^3 



■ (28) 



What we have achieved so far is a classical volume hologram. Given k = 1, this result yields 
the initial quantized signal field a^*^*"^((f) restored with proper amplitude, but quantum 
fluctuations of initial amplitudes of the atoms and the read-out light degrade the fidelity of 
memory below classical limit. The problem comes from the fact that the light and matter 
subsystems do not "forget" their initial states during the interaction. 



III. DOUBLE PASS QUANTUM VOLUME HOLOGRAM 

The double pass operation of our volume hologram looks similar to the case of the QND 
based thin hologram. In order to compensate initial noise terms and to go beyond the 
classical limit, we introduce at the writing stage the second pass through the atomic sample of 
both (the signal and the classical) light waves in the same directions kg and kc- After the first 
pass one has to rotate the spins around the a;-axis by 7r/2 pulse of auxiliary magnetic field, 
and to apply phase shift of 7r/2 to the output signal wave relative the classical driving wave 
within the interaction volume. The transformation of light and matter variables between 
the first and the second pass, which is described by the same equations as before, is given 
by 

a^(^-)(g) = za^(^'°-)(g), x^^^-H^) = Pr^^^g) = 0;^^^'°"*^?). (29) 
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After all three steps of the write-in procedure we arrive to the following in-out relations for 
the double-pass operation: 



W{out) / 



a 
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Xq 



(g) = i(l-/i')a^('")(g)+m( 1- 
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Po 



W(in) 



{q)- 



W{in) 
P2 



{q. 



(31) 

(32) 
(33) 



The double pass read-out procedure looks similar and is described by the Eqs. (30-33), 
where one has to substitute the label W to R. 

We have omitted expressions for higher-order spin amplitudes because the output am- 
plitude of the signal wave by the read-out a^*^°"*)(g) is decoupled from these observables, as 



seen from the Eq. (30). The result (30) also shows that both by the double pass write-in 



and the read-out, the output signal wave amplitude is completely independent of the input 
wave if the coupling constant is chosen to be k = 1, which is one of conditions for efficient 
state exchange between light and matter. In this sense our present model works as the thin 
quantum hologram of p], where by the first pass the input signal amplitude is imprinted 



into the collective spin wave. After the transformation given by the Eq. (29) and the second 



pass, the relevant collective spin amplitude effectively cancels the input signal contribution 
to the output quantized field. 

At the same time, the initial state of light is stored in the collective spin amplitudes 
^w{out) ^w(out) ^ ^ result of the atom- light interaction during two passes with prop- 
erly adjusted coupling constant we have the write stage completed. Atoms and light have 
exchanged their initial quantum states. 

In order to read out the memory one has to repeat the procedure. The output quantized 
field amplitude for the whole write-in and read-out cycle is found by using the W(out) 
amplitudes as the R(in) variables. Finally we arrive to 

■ 3k^ 

Y ' 3 



) = k\2 - k^)a^^'^\q) + i{l - /i2)a«(^")(g-') - ^ 1 - + 



W(in) 

Po 



{q)+ 
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Z-3 

1^ W{in)j 



(34) 



One optimizes the memory performance by choosing the couphng constant k = 1. The 
input and the output variables of the total write-readout protocol of the volume quantum 
hologram are related as 



a 



R{out) I 



(35) 



where f{q) is the added noise term. This transformation is analogous to the one describ- 
ing quantum holographic teleportation of an optical image [151 HB] and the thin quantum 
hologram of [2j . The noise contributions specific for our model of memory is given by 



(36) 



where we assume k = 1. We perform the backward Fourier transform q — )■ p and arrive at 

1 



/l«(-*)(p ) = A^(-)(p ) + Fip), F{p) ^ -^[Fxip) + iFp{p)l 



(37) 



where we apply the same definition of real (Hermitian) noise quadrature amplitudes Fx,p{p) 
as in |2]. Consider orthogonal pixellized spatial modes associated with the field amplitudes 



averaged over the surface Sj of square pixel j 



N of area 5*. The averaged noise 



quadrature amplitudes and the noise covariance matrix are defined via 

1 



x,p 



dpFx,p{p), 



C^(z,j) = (F^(z)F^(j)), C^(^,j) = {Fpit)Fpij)). 
The collective spin noise quadrature amplitudes look like 



-^FxU) = -^ImXr(*"^(j) + -ReP^^'''\j) - - ^ 



ReP. 



W{in) 



V2 



Fp{j) 



o 



W(in). 



6^5 



ImP 



Wi-in) 



(38) 
(39) 

(40) 
(41) 



As seen from the Eq. (14), the real and the imaginary part of the amplitudes X„ can be 
related to real quadrature amplitudes X„ c and Xn,s of the longitudinal collective spin modes 
whose spatial profile is given by 



Z] COS 



2kn 



z] sm 



2k. 
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We arrive to 

ReX^ij) = ^^n,c(j), ImX„(j) = --Lx„,,(j), (42) 

and similar definitions are applied to the amplitudes P„. In the limit of large number of the 
interference layers disused above, \AkzL\ ^ 1, the "sin" and "cos" modes are orthogonal 
and correspond to independent degrees of freedom of the spin coherence. 

The quality of quantum state transfer |?/'(*"')) — > |^(°"*)) is quantified via the fidelity 
parameter F = Assume the input signal field to be in a spatially multimode 

coherent state. For the image field in a coherent state, decomposed over N pixellized modes, 
the fidelity is given by 



1 -1/2 

det{(5,,+C^(^,j)}det{5,,+C^(z,j)} • (43) 



We evaluate the noise covariance matrix assuming all longitudinal and transverse modes of 
the collective atomic spin to be initially in vacuum state, ((^X^*-*"^(j)j ) = ((-PrM*"*"'*(j)) ) = 
1/2, where i = c,s. This yields, 

11. „ /60^^ 



As shown in jTHj, the fidelity of quantum state transfer for simple multipixel arrays scales 
approximately as the A^-th power of the quantity which is called average fidelity per pixel, 
Fav = {FnY^^ ■ For the input signal prepared in a multipixel coherent state we obtain 
the average fidelity per pixel Fav ~ 0.845 for the whole write-in and read-out cycle of our 
quantum hologram. This is significantly higher than the classical benchmark Fav = 1/2 for 
the holographic storage of an optical image in multipixel coherent state, and than the limit 
of Fav = 2/3 on quantum cloning. 

Theoretically, in order to achieve better fidelity one could perform squeezing of each of 



the added noise quadrature amplitudes in (36). Given perfect squeezing, which we do not 



consider here, this would lead to maximum fidelity Fav 1- 

The optimal value of the coupling constant k = 1 corresponds to k = \/2 for the quantum 
volume hologram of [T], which gives there the storage and retrieval efficiency of 0.10 — 0.15. 
This is significantly less than the non-zero quantum capacity limit on the efficiency of 0.5 
and leads, by such a low value of the coupling constant, to performance on a par with 
a classical hologram only. It should be noted that the experimentally achievable atomic 
samples have spatially dependent distribution of atomic density. The required homogeneous 
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coupling constant could be obtained by choosing a proper transversal profile of the classical 
driving wave. 

The resolving power in space of the double-pass quantum hologram can be characterized 
by a number of transverse modes (pixels), effectively stored and retrieved in the memory. 
Similar to the quantum volume hologram of [Ij and the spatially resolving Raman memories 
[TT] . by the write-in and read-out in the same direction the diffraction does not modify the 
light-matter interaction and does not impose limitations on the resolving power which is 
finally determined by the sample geometry: the effectively stored orthogonal modes should 
propagate within the atomic ensemble. Note that the diffraction by propagation at the cell 



length can be compensated by thin lenses, and our definitions (11, 19) of the input and 
output field amplitudes account for this compensation in explicit form. The estimate of the 
number of stored modes given in [1] applies to the scheme considered here as well. 



IV. CONCLUSION 



We have elaborated a novel version of spatially multimode quantum memory for light - 
double pass quantum hologram. Similar to our previous proposals of thin quantum hologram 
[2] and volume quantum hologram [IJ, the present scheme can be considered as an extension 
of classical holograms into the quantum domain. In our scheme, different spatial modes of 
the incoming field are stored in the corresponding orthogonal spatial modes of the long-lived 
collective spin of an ensemble of ultra-cold atoms in absence of spin rotation in external 
magnetic field. The double pass quantum volume hologram inherits some features of the 
thin and the volume quantum hologram. It is able to store transverse modes of an input 
light signal with the same high density as the quantum volume hologram and the Raman 
memories do, and it requires a fixed and relatively low optical depth as the thin hologram 
does. On the other hand, both the write-in and the read-out cycles of the double pass 
quantum hologram require two passes of light, and for a high fidelity storage (with the 
average fidelity per pixel exceeding 0.845), preparation of the collective spin in a quadrature 
squeezed state is also required. Although we considered spin 1/2 atoms, our analysis can be 
in principle applied to alkali atoms provided that the optical detuning significantly exceeds 
the excited state hyperfine sphtting [3]. 
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